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SUMMARY 


The equation of vertical motion for an airplane flying in gusty air 
is simplified in order that its solution is a function of only two 
parameters, namely, the mass parameter of the airplane and the shape of 
the gust the airplane is penetrating. The solutions of the equation are 
presented in the form of charts that can "be used for estimating rapidly 
and easily the acceleration ratios encountered by airplanes with 
different mass parameters penetrating a sharp-edge gust, a gust of 
arbitrary shape, or a triangular gust. 


INTRODUCTION 


For many problems of gust loads it has been found practical to 
calculate the loads as if they were for a rigid airplane restrained in 
pitch under the action of the gust. Although previous investigations 
have presented solutions for particular mass parameters and gust shapes, 
they have not been found sufficiently accurate. In the present paper 
the equation of motion has therefore been solved for a large range of 
mass parameters to obtain charts of airplane reactions to specific 
gusts and to provide a means for obtaining loads on airplane wings for 
arbitrary gust shapes. The charts presented are based on two- 
dimensional unsteady— lift curves and a numerical solution of the 
equation of motion. 


SYMBOLS 



slope of wing lift curve, per radian 


P 


mass density of air, slugs per cubic foot 


Y 


forward velocity, feet per second 


2 
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w 

airplane weight, pounds 

s 

area of wing, square feet 

C 

mean aerodynamic chord, feet 

s 

acceleration due to gravity, feet per second per 
second 

S 

penetration Into gust, chords 

As 

increment of s used in numerical solution 


normalized unsteady— lift function due to penetration 
of a sharp-edge gust 

F 

forcing function in recurrence formula due to 
pentration of gust 

C Iu 

normalized unsteady-lift function for a unit change 
of angle of attack 

V 

transformed unsteady— lift function (l — Cp^s)) 

U 

gust velocity, feet per second 

2 2 > • < 

, . gradient distances of gusts , chords 

^g 

mass parameter ( — \ 

An(s) 

load— factor increment 

^3 

load— factor increment as computed ty the gust— load 
f ormula \ 1 

V 2W / 

An(s)/An s 

acceleration ratio 


acceleration ratio 
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acceleration ratio for penetration of a sharp— edge 
gust 

I integral of acceleration ratio An/An s 

Subscripts: 

max maximum value 



m 


1 


integer denoting m increments of s in recurrence 
formula (s/As) 

denotes point at which solution applies 


DERIVATION AND USE OF CHARTS 


The following assumptions are made for stating the equation of 
motion: 

(1) The gust velocity is uniform across the span and parallel to 
the vertical axis of the airplane at any instant. 

(2) The airplane is in steady level flight prior to entry into 
the gust. 

( 3 ) The airplane can rise but does not pitch under the action of 
the gust . 

(4) The airplane is rigid. 

( 5 ) The lift increment of the horizontal tail due to the gust and 
the airplane motions is negligible as compared to the wing lift 
increment . 


The following equation for describing the vertical motion of the 
airplane is written in terms of acceleration ratios: 


An(s-L) 

^s 



a) 


Mb) ) 

\ U ^cJ d _ _1_ f Sl c T (si - s) ds 

as -Wo ^ 1 


( 1 ) 


The first term on the right represents the force due to the gust and 
the second, the alleviation due to the vertical motion of the airplane. 
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The values of the load— factor Increment £n s and mass param- 
eter Hg may he obtained from the characteristics of the airplane and 
its flight conditions by the following equations: 


and 


^s = 


dC T 


max 


2W 


( 2 ) 


^g = 


2W 


dC T 

p dV ^ 8 


(3) 


Sharp— Edge Gust 

For a sharp-edge gust (fig. 1(a)), equation (l) can be rewritten 


f^il) =c L (s)-i. f 1 c, ( 3l - s) Ml) 
\ /seg s ^g Jo ^ 1 £n s 


ds 


W 


Equation (4) has been solved by a numerical method according to 
reference 1 and the solution is written in the form of a recurrence 
formula. Utilizing the procedure in this reference and making sub- 
stitutions for the unsteady— lift functions gives the numerical solution 
of equation (4) in recurrence form; thus, 



( 5 ) 
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where 


K = c lJ m ) = C L (s) 


m 


g 


I m _ I mr-1 + 


9[*rM + 19 ( — ' ) - 5 ( + ( 


An 


m 


VAns/m-l V ^ 8^2 Wa, 


(6) 


m-3 


Kn = 


\ 2 * 3 , 


(7a) 


Kg =-^-K 1 


g 


(7b) 


k 3 = 


K^ As / 4A^ 2.9 

H g \ 3 24/ 


(7c) 


Ki. 


Kj_ As / ^ 2Ag\ 


124 3 


g 


(7d) 


*5 = 


■^1 As / 4-A-^ 2_ N 

v a V 3 24J 


( 7 ©) 


*6 - 


2 AsK-|_A^ 


3h 


g 


(7f) 
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!*■ AaK]_A^ 



(7g) 


and. 


Kg « 


AaK-)A^ 


(7h) 


A„ = 1 — C^a) 


in 


(8) 


Equation ( 5 ) is evaluated for a range of mass parameters from 10 
to 100 and the results are shown in figure 2. In the calculations , the 
unsteady— lift functions derived for the two-dimensional wing are 
assumed to he made applicable to the finite wing by replacing the slope 
of the lift curve for the two-dimensional wing by the slope of the lift 
curve for the finite wing. These unsteady— lift functions due to 



angle— of— attack change 




for a two-dimensional wing were normalized 


from those given in reference 2 . 


Although equation (4) can be solved in closed form by means of 
operators, the solution is unwieldy and the numerical method was 
adopted for ease in computation. 


The inaccuracies incurred in the solution by the numerical method 
are principally a function of the mass parameter d and the 

o 

increment As. For high values of and small values of As the 

o 

accuracy of the calculation is increased, conversely for low values 

of p and large values of As the accuracy is decreased. When the 
8 

recurrence equation was evaluated, the value of As was appropriately 
changed for the range of p g considered. 
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Arbitrary Gust Shape 


In order to calculate the response to an arbitrary gust shape , the 
response to a sharp^sdge gust may be used as an indie ial response with 
the gust shape as a forcing function in Duhamel f s integral. (See 
reference 3.) Thus, 


£n( si) 

An a 



( S 1 


s) 



(9) 


Various methods are given in reference 4 for evaluating equation (9) 
numerically; however, if the forcing function consists of a series of 
straight lines, equation (9) can he modified by the following procedure; 
For a linear disturbance starting from s = 0 with a slope l/H and 
extending to infinity (fig. 1 (b)), equation (9) can be rewritten as 
follows : 


An( s) 


An = 



(Mil) da 
/ seg 


( 10 ) 


This equation can be evaluated for different mass parameters by merely 

%i(sV 

■ , ^ n 3 / seg 

l/H. In figure 3 the 


integrating the function 


for different values of mass param— 


responses are given for a range of mas 3 


eter and multiplying the result by the constant 

integrations of the j j 

\An s / seg 

parameters. When the ordinates for a curve are multiplied by the 
slope l/H, the resulting response is the solution of equation (10). By 
superposition of a series of functions of the type given by equation (10), 
this approach may be extended to predict the responses to disturbances of 
the type shown in figure l(c) where the disturbance is composed of 
straight-line segments. Figure 3 is used to evaluate this response by 
the following steps: 
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(a) Evaluate the following Integrals In a manner similar to that 
used for equation (10), with the aid of figure 3: 


(b) Superpose the time histories obtained from step (a) so that the 
disturbance function shown in figure l(c) will be formed by the addition 
of the slopes l/H-^ l/Hg, l/H^, and l/H^. Note that in the course 

of this superposition, the functions corresponding to the slopes 
l/Hg, l/H-^, and l/H^ must be displaced by the intervals Hp for the 

slope 1/H 2 , + H 2 for the slope l/H^ , and + H 2 + for the 

slope l/H^. The final response for the disturbance shown in figure 1(c) 
may then be obtained from the following equation: 



— I(s) 
*2 


(ID 


—I(s) 

H 0 



y 


Ms) 

^n s 


= ^-I(a) - ^-I(s - Hp) - ^-l(s - Hp) 



+ - % - a, - h 3 ) - ii(s _ Hl - ^ - h 3 - h 4 ) 


(12) 
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Although equation (12) appears very lengthy, the actual 
computations are relatively simple since all the integrals for a 
particular mass parameter have already been determined (fig- 3 )- The 
values of each integral multiplied by its appropriate constant can be 
tabulated in separate columns displaced by its corresponding interval. 
The addition of the elements in each row of these columns will determine 


the final response 


An(s) 

^3 ' 


Triangular Gust 


Since the triangular gust shown in figure l(d) was used in many 
calculations, a chart was prepared for this specific shape. The 

response to this type of disturbance for a range of mass parameters p 

8 

and gradient distances H was obtained by the method previously 
described for gust shapes approximated by straight lines. Plotting the 

''An 


maximum values of these responses 


An. 


i^nax 


as a function of H with 


M-g as a parameter gives the chart shown in figure 4 . As a matter of 


interest the values of 


Y^s /max 


are plotted as a function of p 


8 


with the gradient distance H as the parameter in figure 5 . When the 
mass parameter of the airplane (equation (3)) and the gradient distance 

H are known, the value of / is readily obtained directly. 

V^s/max 

Evaluating An g by equation ( 2 ) permits the value of An tt[HT in a 


triangular gust to be determined from the acceleration ratio 


V^s/max 


CONCLUDING REMARKS 


The simplified equation of vertical motion of an airplane has been 
solved to obtain charts of acceleration ratio for different gust shapes. 


10 


NACA TN 2036 


The results permit the rapid estimation of maximum acceleration ratios 
for gusts of arbitrary shape. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va. , August 1, 1949 
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(a) Sharp -edge gush. 




U 



(d) Tnongu/ar gush. 


Figure 1-- Gust shapes- 
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Figure 2.- Acceleration ratio for a wing penetrating a sharp-edge gust for various mass parameters. 


NACA TN 2036 












Figure 3-" Integration of acceleration ratio for a ving penetrating a sharp-edge gust for various 

mass parameters. 
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Figure 5* ~ Maximum acceleration ratio for a triangular gust as a function of mass parameter. 





VJl 


NACA TN 2036 









